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Abstract 

We study gluon propagator in Landau gauge with lattice QCD, where we use 
an improved lattice action. The calculation of gluon propagator is performed 
on lattices with the lattice spacing from 0.40fm to 0.24fm and with the lattice 
volume from (2.40fm) 4 to (4.0fm) 4 . We try to fit our results by two different 
ways, in the first one we interpret the calculated gluon propagators as a 
function of the continuum momentum, while in the second we interpret the 
propagators as a function of the lattice momentum. In the both we use models 
which are the same in continuum limit. A qualitative agreement between two 
fittings is found. 

PACS numbers: 11.15Ha, 12.38. Gc, 12.38.Aw, 14.70.Dj 
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1. Introduction 



Lattice QCD enables us to study the nonperturbative nature of QCD from the first 
principle. As space-time is replaced with a discreted lattice, physical results from Monte- 
Carlo simulations can be obtained if the effect of the finite lattice spacing a is under control. 
For this purpose simulations of lattice QCD are performed on large and fine lattices to 
ensure that the effect is small enough. This is expensive. Using an improved action to 
simulate lattice QCD on small and coarse lattices, it is possible to obtain physical results 
with less expense. The latter was proposed long time ago [|l|. It was suggested that one 
can use perturbative theory to improve a lattice action and the effect of the finite lattice 
spacing a can systematically be removed. However, for lattice QCD perturbative series 
converge so slowly that the improvement is not significant. It was pointed out recently 
that perturbative theory can be made more effective by employing the tadpole improvement 
0. With an improved action implemented with the tadepole improvement for lattice QCD 
some physical results are already obtained from simulations on a coarse lattice. For example, 
hadron masses are well determined . 

In this paper we study gluon propagator in Landau gauge with an improved action on 
coarse lattices, where quarks are neglected. Perturbatively a gluon is expected to behave 
like a massless particle. Because of the color confinement of QCD such expectation can not 
be correct in a real world. A nonperturbative study of gluon propagator is required and can 



help us to understand the color confinement. Previous studies [|12| [Tql are undertaken with 
the Wilson action, i.e., an unimproved action, they already show that gluon propagator is 
more complicated than that from perturbative theory. Especially, in Landau gauge it is well 



known that gluon propagator is infrared finite due to Gribov's copies [10], this is also proven 



to be true on lattice 11 



We study gluon propagator on a series of lattices, of which the lattice spacing a is from 
0.40fm and to 0.22fm. At a = 0.40fm we also investigate gluon propagator with lattice 
volume from (2.40fm) 4 to (4.0fm) 4 . This allows us to study the effect of the finite lattice 
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spacing and of the finite volume. The whole calculation is performed with a PC running at 
400MHz. 

In continuum the gluon propagator D°^ v in Landau gauge is characterized by a single 
function D(q 2 ): 

DfM) = (V - ^^Dtf). (1) 

The index a and b is the color index. We define a lattice version of the function D(q 2 ), which 
is denoted as D^, and calculate it. We find that the function Dl can be well interpreted 
as a function of q 2 for q 2 up to3.5(GeV) 2 . We fit our results with a model proposed in 



| l!q| . As our lattice spacing is rather larger, the effect of the finite lattice spacing needs to be 
studied in detail. For this we also try to fit our propagator as a function of a modified lattice 
momentum variable, where the same model is used with a slight modification suggested by 
perturbative theory. An qualitative agreement is found between these two fittings. 

Our paper is organized as the following: In Sect. 2 we introduce the action and our 
notation, we discuss different interpretations of our data. In Sect. 3 we fit our results by 
taking D/, as a function of q 2 with the model in fll2"]. In Sect. 4 we fit our results by taking 



-Dl as a function of lattice momentum with the same model, but with a slight modification 
as suggested by perturbative theory. In the both sections we will discuss the effect of the 
finite lattice spacing and of the finite volume on our fitting results. Sect. 5 is our conclusion. 

2. The Action and Our Notations 

We take the one-loop improved action for gluon @ , where the action consists of plaquette, 
rectangle and paralellogram terms and is accurate up to errors of 0(a 2 a 2 , a 4 ). Implementing 
tadpole improvement the action becomes || 

S(U) =/?E ~ReTr(l - U pl ) + (3 rt £ iReTr(l - U tr ) + (3 pg ]T ^ReTr(l - U pg ), 

P l 6 rt 6 pg 6 

Prt = -^(1 + 0.4805a s ), /3 M = -4o.03325a s , 

M = (-ReTr(tfpi))*, a s = ^ ^—^L 2 

v 3 \ pi/j , 3.06839 V ; 
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We used this action to generate gluonic configurations at j3 = 6.8, (3 = 7.1 and (3 = 7.4, 
where the lattice size is 6 4 , 8 4 and 10 4 respectively. The parameter w is determined by self- 
consistency. We also generate configurations at f3 = 6.8 with lattice size of 8 4 and 10 4 . The 
configurations are generated with the pseudo heat bath method ||, where the three SU(3) 
subgroups were updated 3 times in each overall update step. Configurations are separated 
by 40 sweeps to ensure that they are statistically independent. The configurations used in 
this work are summarized in Table 1. 



Table 1 



p 


L 4 


No. 


a 


Uq 


6.8 


6 4 


500 


0.40fm 


0.8267 


6.8 


8 4 


200 


0.40fm 


0.8267 


6.8 


10 4 


200 


0.40fm 


0.8267 


7.1 


8 4 


500 


0.33fm 


0.8434 


7.4 


10 4 


250 


0.24fm 


0.8631 



In Table 1 the lattice spacing a is determined by the string tension ||. We use the naive 
steepest descent methode to fix the gauge of our configurations 0. We define: 

A(x) = '^2\U u (x — av) — U u (x) — h.c — trace term], 

9 = ^H Tr ( A ( x ) At ( x ))> (3) 

where V = L 4 is the lattice volume. By gauge transformation we try to minimize the 
quantity 6.9 = corresponds to the Landau gauge. Numerically we require 6 < 10~ 8 . The 
gauge field is defined as 

A,(x + aii/2) = ^(f/,(x) - Ul(x)) - JL-Ttp^x) - U^x)). (4) 

It should be noted that in the Wilson action the bar coupling g is related to the (3 as 
g P = 6. With the action in Eq.(2) the bar coupling g is related to (3 as 
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% = j, (5) 

if one neglects the tadpole improvement and the correction in the action at one-loop level. 
Adding the improvement and the correction the relation becomes complicated. In Landau 
gauge the quantity Z) x x) is independent on t for each configuration. We checked this 
under the numerical condition for 6. The quantity varies at different t under 0.1%. 
The Fourier transformed gauge field is 

A,(q)=a 4 Y,e- ig - {x+afl/2) A,(x), (6) 

X 

where q is the momentum: 

% = -j^n^ n M = 0,l,---,L-l. (7) 

We define a function on lattice as 

D L(Q) = < 0|Tr(^(g)4(g))|0 > . (8) 

In the limit of a — > the function Di(q) will approach to the function D(q 2 ) in the 
sence that it will be proportional to D(q 2 ). We will calculate the function D L (q) with the 
configurations listed in Table 1. For further convenience we introduce 

2 1 

g„ = - sin(-ag M ), q 2 = J2 ?2> ^ = 2 «J- ( 9 ) 
a z ^ ^ 

In general the calculated function D L (q) will depend on g 2 , g 4 and other possible in- 
variants on lattice, if ag M is small enough, Di,(q) will approximately depend only on the 
variable 

Previous calculation on large lattice shows that even for a 2 q 2 < 1 D L (q) still can not be 
interpreted as a function of q 2 , the data points for D^q) with momenta directed along 
one of the four axes behave differently than those with momenta directed off axis and the 
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difference is significant. This clearly indicates that the rotation invariance is violated due to 
finite lattice spacing and the effect is large. With a improved action like the one in Eq.(2) 
the effect of the violation will be reduced. Our data for D^q) from the configurations with 
(3 = 7.4 is shown in Fig. la, where Dl{q) is plotted as a function of q 2 , and data with same 
q 2 are averaged. From the figure we see that the data points behave well as a function 
of q 2 in a range from a 2 q 2 = to a 2 q 2 = 5.13. For a 2 q 2 > 5.13 the effect due to the 
violation of the rotation invariance becomes significant, indicated by a small jump in our 
data arounda 2 g 2 = 6.3. This is clearly shown by the data from our smallst lattice in Fig. lb, 
where a small jump is roughly at a 2 q 2 = 10.0 or at q 2 = 2.5GeV 2 , for a 2 q 2 > 8.0 the data 
points become unregular. The similar behave is also found for the lattice with the size L = 8. 
We plot our data from the lattice with (3 = 7.1 in Fig.lc. We will try to fit our data lying 
between q 2 = and the q 2 , where the data points begin to become unregular, by taking Dl 
as a function of q 2 . 

However it may be wondered that the effect of the finite lattice spacing is small in the 
above consideration. To study this we also try to fit our results by taking D L (q) as a function 
of q 2 . The replacement of q 2 by q 2 can be thought as kinematic correction. In Fig. 2a we plot 
Di(q) from the lattice with f3 = 7.4 as a function of q 2 , where data points with adjacent 
momenta lying within 5a 2 q 2 < 0.0001 are averaged. From Fig. 2a one can see that the data is 
not well scaled by q 2 and this indicates that Dl may depends on on other lattice invariants 
and the kinematic correction needs to be modified. Perturbative calculation with the action 
in Eq.(2) also shows that D^{q) does not only depend on q 2 but also on q A . At the tree-level 
without the tadpole improvement one can obtain [I7j 



where there is a constant term. With this in mind the kinematic correction above needs to 
be further modified. We introduce a new variable: 

?L = ? + ^<f, (12) 
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and for a — > 0: 



gl = g 2 + 0(a 4 ). (13) 

In Fig. 2b we plot the gluon propagator as a function of q\, where data points with 
adjacent momenta lying within 5a 2 q 2 L < 0.0001 are averaged. From the figure one can see 
that Dl can be interpreted as a function of q 2 L) although there is still a small fluctuation 
among the data points, but it is much better than interpreting Dl as a function of q 2 . We 
will also fit our results by taking Dl as a function of q L . 

It should be pointed out that only on-shell quantities calculated with an improved action 
are improved. For an operator in general one needs to improve it to remove effect of finite 
a. In our case it is complicated to do so. The reason is that the primary field on lattice is 
the gauge link U^ix), certain effect of the finite a at order of a 2 is introduced by extracting 
the gauge field A^(x) from U^(x), and some effect at order of a 2 arises if one calculate 
the propagator. To obtain an improved -Dl(qO one should improve the extraction and then 
improve the operator A tl (q)A^ (q). In this work we will not consider to obtain an improved 
D L (q). 



3. Fitting propagators as a function of continuum momentum 

In this section we fit our calculated propagators as functions of q 2 with the model |L2| 



D ^ 2 ) = {M2 yJ + (g2)1+a (14) 

where M is a dimensional parameter and a is the anomalous dimension. Z is a parameter 
related to go and the renormalization constant of wave-function. In this and the next section 
we will not discuss this parameter. In the fitting of our propagators we exclude the data 
points with the lowestg 2 including the point with q 2 = 0. We try to include more data points 
with larger q 2 untile \ 2 P er degree of freedom close to 1, but not larger than 1. Beside this 
the total fitting range is limited because of the reason discussed in the last section. Our 
fitting results are summarized in Table 2: 



Table 2 





L 


Range 


Z 


a 2 M 2 


a 


X 2 /dof 


6.8 


8 


3-9 


5.811(133) 


1.654(92) 


0.107(24) 


0.35 


6.8 


10 


3-12 


5.975(77) 


1.730(49) 


0.111(14) 


0.07 


7.1 


8 


3-12 


5.511(14) 


0.949(106) 


0.187(18) 


0.28 


7.4 


10 


3-14 


5.298(75) 


0.602(41) 


0.312(13) 


0.08 



The data points are labeled by an integer n, smaller n corresponds to smaller q 2 , n = 1 
corresponds to q 2 = 0. The range is the fitting range and given by this integer. We do 
not list our results for the configurations with (3 = 6.8, L = 6, because the fitting quality 
is poor with small number of data points. We also tried to fit our data with other models. 
For example, the model proposed in ||15|| , but it can not well describe our data in the range 



given above by the indication of large x 2 /dof. To illustrate our fitting quality we plot our 
data with fitting results in Fig.3a and in Fig.3b for (3 = 7.4 and for (3 = 6.8 respectively. 

From our fitting results and with the lattice spacing given in Table 1 we obtain the 
dimensional parameter M in GeV: 

M = 0.633(28) for (3 = 6.8, L = 8, M — 0.647(14) for (3 = 6.8, L = 10, 

M = 0.582(56) for (3 — 7.1, L = 8, M — 0.636(34) for (3 = 7.4, L = 10. (15) 

We note that the values of M at /3 = 6.8 with L = 8 and with L = 10 can be taken as the same 
within the errors, this may indicate that the effect of the finite volume is small. However, we 
have only two values from two lattices, where the difference in volume is not large, a study 
on larger lattices with different volumes is needed to confirm the above conclusion. From 
above the values of M remain relatively stable with different lattice spacings a. We use a 
standard weighted least-squares procedure to average the M- values and obtain the central 
value M c : 

M c = 641(ll)MeV. (16) 



We also note that the deviation of M at f3 — 7.1 from M c is the largest, but it is still in 
the range of the errors. The deviation of other values is small. This leads to conclude that 
there is no significant dependence of M on the finite a and to interpret the central value M c 
as the parameter M in the continuum limit. 

As the next we study the anomalous dimension a. From the fitting results in Table. 2, 
the values from two lattices with different volumes at (5 — 6.8 are the same within the errors. 
This may also indicate that the effect of the finite volume in a is small as that in M. For 
values at different (3 there is clearly a a-dependence of a. We try to model the dependence 
with a simple formula: 

a (a) = «o + ct\a 2 (17) 

where a is in unit of fm and a can be thought as the anomalous dimension in continuum. 
At (5 — 6.8 we take the value of a with L — 10 for the fit. The results is 

a = 0.421(20), «! = -1.97(19), x 2 = 1-42. (18) 

In Fig. 4 we illustrate the quality of this fit. The value of x 2 is acceptable. With the above 
discussion we conclude in this section with our data that the gluon propagator in continuum 
has a mass scale at 641(ll)MeV and there is an anomalous dimension which is 0.421(20). 

4. Fitting propagator as a function of lattice momentum 

In this section we fit our data with the same model as in the last section, where the 
variable q 2 is replaced by g|, and a constant term is added as suggested by perturbative 
theory. This term will vanish in the limit of a — > 0. The model with such modification is: 

^)= {Ap)1 J +i - ly+a + ^. (19) 

The data selection for fitting is similar as in the last section. We discarded the first two 
points and try to include more data points in the direction of increasing q 2 , until x 2 P er 
degree of freedom close to 1 but never larger than 1. Our fitting results are summarized in 
Table 3. 



Table 3 



p 


L 


Range 


Z 


a 2 M 2 


a 


c 


X 2 /dof 


6.8 


6 


3-12 


6.297(72) 


1.433(67) 


-0.019(8) 


-0.319(11) 


0.65 


6.8 


8 


3-13 


6.204(68) 


1.496(56) 


-0.080(10) 


-0.379(14) 


0.50 


6.8 


10 


3-16 


6.007(59) 


1.440(46) 


-0.097(11) 


-0.365(14) 


0.44 


7.1 


8 


3-13 


6.203(69) 


0.945(39) 


0.026(12) 


-0.383(15) 


0.76 


7.4 


10 


3-15 


5.809(56) 


0.546(20) 


0.147(13) 


-0.361(16) 


0.58 



The data points are labeled by an integer n, smaller n corresponds to smaller q 2 Ll n = 1 
corresponds to q\ — 0. The range is the fitting range and is given by this integer. To illustrate 
out fitting we plot our data with the fitted curves in Fig.5a and Fig. 5b as examples. From 
above one can realize that the c roughly remain as the same at different (3 and at different 
lattice volumes except the one at f3 — 6.8, L = 6. For other values of c the variation is in 
the range within errors. This fact will ensure that the term with c in Eq.(19) will vanish in 
the limit of a — > 0. 

With the lattice spacing in Table 1 we obtain the mass parameter: 

M = 0.590(11) for (3 — 6.8, L = 6, 

M = 0.602(9) for (3 = 6.8, L = 8, M — 0.591(8) for (3 = 6.8, L = 10, 

M = 0.580(12) for (3 — 7.1, L = 8, M — 0.606(42) for (3 = 7.4, L = 10. (20) 

We observe that the value of M is within the errors stable in the range of (3 and in the 
range of the volume of our lattices. This also indicates as in the last section that the gluon 
propagator in continuum limit contains a mass scale. We calculate with the data the central 
value: 



M c = 592(5)MeV. (21) 
This value is 10% smaller than the one in the last section. The reason may be due to the 
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difference between the variables q 2 and q\, i.e., the effect of the finite lattice spacing is still 
observable, it may be at order of 10%. 

At j3 = 6.8 the parameter a from the lattice with the volume L = 6 is quite different than 
the values from other larger lattices, while there is no large deviation between the values 
from L = 8 and from L = 10. The similar is also found with the parameter c. This may be 
explained that at (3 = 6.8 the effect of the finite volume is significant for L smaller than 8. 
We try to fit this dependence with the same relation as in the last section: 

a = ceo + aid 2 , (22) 

the fitting results are: 

a = 0.285(20), e*i = -2.38(17), \ 2 = 0.07. (23) 

The fitting quality is better than that in the last section. But «o is different than this in 
the last section. Here cto may not be interpreted as the anomalous dimension. The reason 
is: In the model in Eq.(14) the effect of the finite a is only contained in the parameters, 
while in the model in Eq.(19) not only the parameters implicitly contain possible effect of 
the finite a, but also some effect of the finite a is explicitly included, representing by the 
term with c and with the variable q\. If the latter can be neglected, one may interpret «o 
as the anomalous dimension, which is 30% smaller than ao in the last section. This may be 
thought as a qualitative agreement. 



5. Conclusion 

In this work we studied gluon propagator with an improved action on lattice. Unlike 
with the Wilson action, we find that the effect due to the violation of the rotation invariance 
is reduced, resulting in that D L can be well interpreted as a function of q 2 as a continuum 
variable for a 2 q 2 not larger than 5.13. We fitted the gluon propagator with a model proposed 



m 



12| and studied the effect of the finite lattice spacing and of the finite volume. Our 



conclusion is that in the continuum limit the gluon propagator has a mass scale around 
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600MeV and a anomalous dimension which is 0.412. We also tried to study the gluon 
propagators as a function of an improved lattice momentum q\ and fitted them with the 
model which is the same in the continuum limit as the one used before. The fitting results 
are qualitatively in agreement with the study of as a function of q 2 . The difference 
between two fittings may be explained by the effect of the finite lattice spacing, and by the 
effect depending on how to model the gluon propagator Dl- Simulations with large lattices 
may be required to study the difference in detail. 

Finally it should be pointed that calculations with larger lattice than those in our work 
will allow us not only to study the effect of the finite lattice spacing and of the finite volume 
in more detail but also to have more data points with small q 2 . 
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Figure Caption 

Fig. la: The gluon propagator D L at f3 — 7.4 with L — 10 plotted as a function of q 2 . 

Statistical errors are not plotted, they are small and around 1%. 

Fig. lb: The gluon propagator D L at (3 — 6.8 with L = 6 plotted as a function of q 2 . 

Fig.lc: The gluon propagator D L at f3 — 7.1 with L = 8 plotted as a function of q 2 . 

Fig. 2a: The gluon propagator Dl at (5 — 7.4 with L = 8 plotted as a function of q 2 . The 

x-axis is for a 2 q 2 . 

Fig. 2b: The gluon propagator D L at f3 — 7.4 with L = 8 plotted as a function of q\. The 
x-axis is for a 2 q 2 L . 

Fig. 3a: The gluon propagator D L at f3 — 7.4 with L = 10 plotted as a function of q 2 , the 
line is from our fit. All data shown in this figure are used in fitting. 

Fig. 3b: The gluon propagator D L at (5 — 6.8 with L = 8 plotted as a function of q 2 , the line 

is from our fit. All data shown in this figure are used in fitting. 

Fig. 4 The extrapolation of a. The x-axis is for a 2 in (fm) 2 , the y-axis is for a. 

Fig. 5a: The gluon propagator D L at (5 — 7.4 with L = 10 plotted as a function of q 2 Ll the 

line is from our fit with the model in Eq.(19). All data shown in this figure are used in 

fitting. 

Fig. 5b: The gluon propagator D L at (3 — 6.8 with L = 6 plotted as a function of q 2 L) the line 
is from our fit with the model in Eq.(19). All data shown in this figure are used in fitting. 
Fig. 6 The extrapolation of a. The x-axis is for a 2 in (fm) 2 , the y-axis is for a. The a is 
defined in Eq.(19). 
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